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A search for a sidereal modulation in the MINOS near detector neutrino data was performed. If 
present, this signature could be a consequence of Lorentz and CPT violation as predicted by a class 
of extensions to the Standard Model. No evidence for a sidereal signal in the data set was found. 
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implying that there is no significant change in neutrino propagation that depends on the direction of 
the neutrino beam in a sun-centered inertia! frame. Upper limits on the magnitudes of the Lorentz 
and CPT violating terms in these extensions to the Standard Model lie between 0.01-1% of the 
maximum expected, assuming a suppression of these signatures by factor of 10~^^. 

PACS numbers: 11.30.Cp,14.60.Pq 



At experimentally accessible energies, signals for 
Lorentz and CPT violation can be described by a class 
of extensions to the Standard Model often referred to 
as the SME [l], This SME theory is an observer- 
independent framework that contains all the Lorentz vio- 
lating (LV) terms involving particle fields in the Standard 
Model of particle physics and gravitational fields in the 
General Theory of Relativity. Since the Standard Model 
is thought to be the low-energy limit of a more funda- 
mental theory that unifies quantum physics and gravity 
at the Planck scale, mp ~ 10"'^^ GeV, it is suggested 
in that the violations of Lorentz and CPT symme- 
tries predicted by the SME provide a link to Planck scale 
physics. Although the magnitude of LV signatures in the 
accessible energy limit are suppressed by a factor of or- 
der the electroweak scale divided by the Planck scale, 
mw/mp - 10-1^ [3, these low-energy probes of new 
physics have been explored in many ways with current 
experimental technologies 

The SME framework predicts several unconventional 
phenomena, among which is one that arises from the de- 
pendence of the neutrino oscillation probability on the 
direction of neutrino propagation [1, . For experiments 
like MINOS @] with both beam neutrino source and de- 
tector fixed on the Earth's surface, the Earth's sidereal 
rotation causes the direction of neutrino propagation p to 
change with respect to the Sun-centered inertial frame in 
which the SME is formulated 0]. The theory predicts 
that this rotation introduces a sidereal variation in the 
number of neutrinos detected from the beam. The LSND 
collaboration Q did not see this signal. In this paper we 
use a sample of neutrinos identified in the MINOS near 
detector (ND) in a search for this sidereal signal. The 
neutrinos were generated by the Neutrinos at the Main 
Injector (NuMI) neutrino beam at Fermilab Q, a beam 
whose flavor composition is 98.7% Vf^-V v^^ 

According to the SME, the probability that one of 
these oscillates to flavor u^^ where a; is e or r, over 
a distance L from its production to detection due to 
Lorentz and CPT violation is given by Q 

~ L\C)^:^ + {Ac)^x cos (weTe) 

+ {As)^tx sin {lO(^T(^) + {Bc)fix cos (2a;0T®) 
+ (Z?,)^,sin(2c^e^e)]^ (1) 

where cj® = 27r/(23'' 56" 04.0982") is the Earth's side- 
real frequency and T® is the local sidereal time of the 
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neutrino detection. For the MINOS ND, < L 750 m. 
In this equation, the expressions for {Ac)^_ix, {As)^ix, and 
{C)^x include both CPT and Lorentz violating terms; the 
expressions for {Bc)^x and {Bs)p,x include only Lorentz 
violating terms. Since CPT violation implies Lorentz vi- 
olation in field theory ^ , there are no terms that depend 
on CPT violation alone. These parameters are combi- 
nations of the SME coefficients (a^)^ and {clY'^ that 
describe LV Q. The magnitude of the coefficients in 
eq. ([T]) also depend on the neutrino propagation direc- 
tion. For the NuMI beam line, the direction vectors of 
the beam are defined by colatitude x = (90°— latitude) 
= 90.° - 41.84056333°, 9 = 93.346°, and = 23.909°. 
Here 9 and (j) are the beam zenith and azimuthal angles; 
9 is measured from the z-axis which points toward the 
zenith; </> is defined from the x-axis which is along the 
detector axis and increases away from the NuMI target. 
The y-axis makes a right handed coordinate system. The 
explicit relationships between the SME coefficients and 
the beam direction are found in fsj. In eq.([T|), the CPT 
violating terms depend on L and the Lorentz violating 
terms depend on L x where Ey is the neutrino en- 
ergy. This unconventional behavior is to be compared 
with the LjEv dependence of the oscillation probability 
resulting from non-zero neutrino masses p!o| . 

The MINOS ND pjj] is a magnetized 0.98 kton steel 
scintillator tracking calorimeter that lies 103 m under- 
ground at Fermilab and it is made of 282 4 x 6 vc? oc- 
tagonal planes. Each plane is comprised of a 2.54 cm 
thick steel plate and a layer of scintillator strips with di- 
mensions of 4.1 cm X 1 cm. Each scintillator strip is 
coupled via wavelength-shifting fiber to one pixel of a 
64-pixel Hamamatsu M64 PMT [T5|. The ND readout 
continuously integrates the PMT charges, while the data 
acquisition accepts data above a threshold of 0.25 photo- 
electrons. 

Protons with mean energy of 120 GeV/c are extracted 
from the Main Injector giving a beam spill of 10 /zs du- 
ration and period of approximately 2.4 s. The number 
of protons delivered to the target for each spill (POT) 
was measured using toroidal beam current transform- 
ers. The uncertainty in the number of POT for each 
spill is ±1.0% [6j]. Neutrino events in the ND fiducial 
volume were selected based on their timing and spatial 
information [l^ . The neutrinos were separated into 
charged current-like (CC) or neutral current-like (NC) 
events, as described in Q. CC events, identified by a 
/Lt~ track were selected for further analysis to maxi- 
mize the — » Vx oscillation signal in MINOS. 

Standard beam quality cuts @ were applied to select 
spills for the analysis. In addition, data quality cuts were 
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TABLE I: Run Parameters 





CC Events 


POT 


Run Dates 


Run I 
Run II 


1.82 X 10" 
1.62 X 10** 


1.25 X 10^" 
1.14 X 10^° 


May05 - FebOG 
Sept06 - Mar07 



applied to remove runs in which there were detector prob- 
lems, including cooling system failures, magnetic coil fail- 
ures, or an incorrectly configured readout trigger. 

The data were taken during two run periods. The pa- 
rameters for these two runs are given in Table [H The 
numbers of events and POT given are the numbers re- 
maining in the sample after all cuts have been made. 

Since the sidereal phase histograms in this analysis re- 
quire accurate event timing, we describe how time stamps 
are generated. The spill time is determined by the Global 
Positioning System (GPS) receiver located in the ND 
hall that reads out absolute Universal Coordinated Time 
(UTC) and is accurate to 200 ns The Main Injec- 

tor accelerates protons to 120 GeV/c and the spills are 
extracted to NuMI using a pulsed dipole magnet. The 
GPS time of the extraction magnet signal is recorded 
and defines the spill time p^ . 

Each neutrino event was tagged with the local sidereal 
time (LST) of its spill - the GPS spill time converted 
to sidereal time. The local sidereal phase of an event is 
given by LSTx(cL'0/27r) and has a range of 0-1. Event 
times were not corrected for their time within a spill, an 
approximation that introduces no significant systematic 
error into the analysis. 

The events in each spill were placed into a single bin in 
a histogram spanning 0-1 in local sidereal phase (LSP). 
The POT in the spill were binned into a second LSP 
histogram. By dividing these two histograms, we get the 
number of i/fj, events/POT as a function of LSP. This 
final histogram gives the normalized neutrino event rate 
in which we search for sidereal variations. 

We used 32 bins for the LSP histograms. This binning 
was chosen to search for sidereal variations with a Fast 
Fourier Transform (FFT) [l^ and the algorithm works 
most efficiently for 2^ bins. Since eq.((T]) puts power into 
harmonic frequencies associated with Fourier terms to 
nw®, and for this analysis n = 1-4, we chose A/" = 5 as 
the minimal binning that retains these harmonic terms. 
Each bin spans 0.031 in LSP or 45 min in sidereal time. 
The histograms of the i/^ events/POT as a function of 
LSP for Run I and Run II are given in Fig. [TJ The dif- 
ferences in the average event rates are due primarily to 
different relative positions of the target and magnetic fo- 
cusing horns for the two runs. 

We performed an FFT analysis on the Run I and Run 
II sidereal phase histograms in Fig. [1] and we computed 
the weighted mean of the powers returned for the even 
(cos ) and odd (sin ) powers for harmonic frequencies out 
to 4a'0T0. The weighting factors were the mean event 
rates for each run. The resulting mean powers, p{FFT), 
are listed in Table |TT1 
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FIG. 1: The local sidereal phase histograms for Run I and 
Run II. Superposed are fits to a constant sidereal rate. 



TABLE II: Weighted mean of Run I and Run II FFT powers 
in first four even/odd harmonic coefficients; Vf is the proba- 
bility that the mean power is a noise fiuctuation. 



cos p(FFT) Vf 


sin p(FFT) Vf 


(aj(Br®) -0.002 0.91 
(2tj®T©) 0.011 0.54 
(3a;®Te) -0.006 0.74 
(4tj®Te) -0.016 0.37 


(l^cbT®) 0.024 0.18 
(2a;®T©) 0.011 0.54 
(3a;®r®) -0.004 0.83 
(4a;®T®) 0.023 0.20 



These results were tested for several possible system- 
atic effects. We found that systematic increases or de- 
creases in the event rate of 5% in 6 months do not affect 
these results. We also searched for systematic changes 
in the rates from day to night and found no variations 
> 0.1%. In addition, we searched and found no sidereal 
modulation in the CC/NC ratio for these data. This test 
shows that there are no systematic effects associated with 
neutrino production in the beam that affect this sidereal 
analysis. 

We constructed 1,000 simulated experiments for both 
runs without a sidereal signal to test the significance of 
the powers given in Table [Til We used the data them- 
selves to construct these experiments. We first generated 
1,000 sets of sidereal phases for Run I and Run II, with 
each set having the same number of entries as spills in the 
run. The phases were drawn randomly from the sidereal 
phase distribution constructed from the start times for 
each spill in the Run I or Run II data set. We then put 
the events for each spill in the run into 1,000 separate 
histograms according to the scrambled sidereal phase as- 
signed. The number of POT for that spill was entered 
into a second set of 1,000 histograms according to the 
same set of sidereal phases. The division of each event 
histogram by its corresponding POT histogram results in 
the simulated experiments - histograms of the number of 
events/POT as a function of LSP without a signal. 
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TABLE III: Limits to SME coefficients for — » in terms 
of the suppression factor mw /rnp ~ 10"^^; have units of 
(GeV) and cl are unitless. 



-0.15 -0.1 -0.05 0.05 0.1 0.15 
FFT Power 



FIG. 2: The distributions for the even (cos ) and the odd 
(sin ) mean powers for harmonic frequencies to 'iuj^T^, from 
the FFT analysis of 1,000 simulated experiments in Run I and 
Run II. Superposed on these distributions is a Gaussian fit of 
width a = 1.8 X 10"^ This fit was obtained independently 
for both distributions. Values outside of the vertical lines are 
more than Sa from the mean. 



We performed the same FFT analysis on each of the 
1,000 Run I and 1,000 Run II simulated experiments as 
was done with the sidereal phase histograms in Fig. [T] 
The powers returned by these FFTs give the fluctua- 
tion spectrum expected from sidereal phase histograms 
in which there is no sidereal signal. As for the data, we 
computed the weighted mean power for each harmonic 
in a pair of simulated Run I and Run II experiments. 
The distributions for the even (cos ) and the odd (sin ) 
mean powers for harmonic frequencies out to 4a;^T0 are 
shown superposed in Fig. [2l Clearly these even and odd 
distributions arc nearly identical. In addition, a Gaus- 
sian of width (T = 1.8 X 10~^ has been superposed onto 
these two distributions in Fig. [51 This fit was obtained 
independently for both distributions. We use this Gaus- 
sian to estimate the probability that the powers returned 
by the FFT analysis of the sidereal phase histograms in 
Fig. [T] are due to statistical fluctuations. 

Tablellllgives the probability, Vp, that the mean power 
represents a noise fluctuation. It was calculated as the 
probability of drawing a value of the weighted mean 
power for the two data sets at least as large as found 
from the parent Gaussian distribution in Fig. [2l Since 
the largest fluctuation in the FFT power in the Fig. [1] 
histogram is 1.32cr we conclude that no term reaches the 
level of a 3cr detection. We have determined that these 
results are insensitive to the exact choice of the zero point 
of sidereal phase. This model-independent result implies 
that there is no signiflcant change in normalized neutrino 
event rate that depends on the direction of the neutrino 
beam in a sun-centered inertial frame. In the context of 
the SME, this result is inconsistent with the detection of 



xlO"" 


xlO"" 


at 3.0 X 10"^ 
cl^ 0.9 X IQ-^ 
cf^ 5.6 X 10"* 
2.7 X 10"* 
^ 1.3 X 10"* 


al 3.0 X 10"^ 
cJJ^ 0.9 X lO"'' 
cV 5.5 X 10-* 
cr 1.2 X 10-* 



LV. 

In the absence of a sidereal signal, we can establish up- 
per limits on the SME coefficients (a^)'' and [clY'^ that 
describe LV Q using the standard MINOS Monte Carlo 
simulation. The simulation includes weightingto account 
for hadron production off the NuMI target [(|. In this 
simulation, events are generated by modeling the NuMI 
beam line, including the hadron production by the 120 
GeV/c protons on target, the propagation of the hadrons 
through the focusing elements and 675 m decay pipe to 
the beam absorber, and the calculation of the probabil- 
ity that any neutrinos generated traverse the ND. The 
ND neutrino event simulation takes the neutrinos from 
the NuMI simulation, along with an energy determined 
by decay kinematics, and uses this information as in- 
put into the simulation of the ND. With the known L 
and El, for the simulated neutrino events, as well as the 
beam direction, we can inject a Lorentz-violating signal 
into cq.(IT]). The construction of MC-generated sidereal 
phase histrograms is described elsewhere [l6l |. 

The limits on the LV coefficients [ah)^ and {clY'"'^ 
were determined from a set of 200 simulated experi- 
ments. First we set all but one LV coefficient to zero. 
We next weighted the simulated neutrino events in each 
histogram by its survival probability computed accord- 
ing to eq.lHI), assuming the LV coefficient is small. We 
then increased the magnitude of the nonzero coefficient 
until one of the FFT powers in the simulated phase his- 
togram was 3(7 away from the mean of the distribution 
in Fig. [21 An average of the 200 determinations of each 
SME coefficient, scaled in terms of the suppression factor 
mvv/iTip 10"^^, is given in Table [Till This procedure, 
by which we vary one parameter at a time to determine 
the limits, could miss fortuitous cancellations of SME co- 
efficients thereby masking a signal. However, we consider 
such cancellations in Nature to be highly unlikely. 

In summary, we find no significant evidence for side- 
real variations in the MINOS ND neutrino data. When 
framed in the SME theory [5], this result leads to the 
conclusion that we have detected no evidence for the vio- 
lation of Lorentz and CPT invariance. Based on these re- 
sults, we computed limits on the LV SME coefficients and 
find that their magnitude is < 1% of the suppression fac- 
tor mw/mp ^ 10^^''. For the a^-type SME coefficients, 
the MINOS limits are a factor of 3 lower than those re- 
ported by LSND f^; for the CL-type SME coefficients, the 
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MINOS limits are at least 4 orders of magnitude lower 
than LSND's. 
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